12.1 FortranとCの混合
Fortran ソースコードとCのコードを混在させる場合でFortran 2003のC相互運用を利用しない場合には、 以下の点に留意してください。
12.1.1 ネーミング規約
外部手続き名と共通ブロック名は小文字にアンダスコアを付加したものです(これは 標準UNIX f77と同じです)。メインプログラム単位の名前は‘main
’と
なります。
モジュール変数とモジュール手続き名はモジュール名(小文字)に
‘_MP_
’を付加、さらに言語要素名(小文字)を付加する形で構成され
ます。
注:これらの規約は-compatible オプションを使用した場合には異なる可能性があります。 Windowsにおいてこのオプションを指定すると、 外部手続き名と共通ブロック名は大文字でかつアンダスコアが付かないものとなる一方、 モジュール要素はモジュール名を小文字ではなく大文字で持つことになります。 (64ビットシステムでは常に-compatibleです。) その他の多くの環境では、 このオプションの使用により名前にアンダースコアが既に含まれている場合に、 (通常付け加えられるもの後に更に)アンダースコアが追加されます。
12.1.2 初期化と終了
メインプログラムがFortranで書かれていない場合には初期化ルーチン__nAGf90_rts_init(int argc,char*argv[])または入出力初期化ルーチン
__nAGf90_io_init(void)のいずれかをCからコールする必要があります。
__nAGf90_rts_init
ルーチンはFortranの浮動小数点環境を初期化すると共に、
Fortran 2003の組込関数及びF90_UNIX
モジュールからのコマンド行引数に対するアクセスを可能にします。
またプログラムは__nAGf90_finish(int status)
によって終了させる必要がありま
すが、通常のC終了ルーチンを用いてすべてのFortranファイルのクローズ、Fortran
入出力バッファのフラッシュを行う前に、__nAGf90_io_finish(void)
をコールする
というアプローチも取れます。(フラッシュを行わないとオープンされている出力ファイルが壊れる可能性があります)
12.1.3 呼出し規約
以下のセクションでは、nAG Fortranコンパイラを用いてコンパイルされたFortran プログラムによって用いられる呼出し規約を、Cの用語を用いて詳細に記述します。 従ってこの情報はFortranとCのコードを混在させたいと考えている人にとって最も 有用と言えます。
近代FortranのFortran 77サブセットを対象とした規約はデファクトである
Unix f77規約と互換性があります(-f77 オプションを用いずに
コンパイルされたCOMPLEX
関数を除く)。
12.1.4 データタイプ
CからFortranへの交信に役立つデータタイプの定義は、通常/usr/local/lib/nAG_Fortran
(UNIX系のシステムの場合)に置かれているdope.h
とnagfortran.h
ファイルの中にあります。
Fortranデータタイプ | Fortran精度 | C typedef 名 |
INTEGER | 8 bits | Integer1 |
16 bits | Integer2 | |
32 bits | Integer または Integer3 | |
64 bits | Integer4 | |
LOGICAL | 8 bits | Logical1 |
16 bits | Logical2 | |
32 bits | Logical | |
64 bits | Logical4 | |
REAL | single | Real |
double | Double | |
double-double | DDReal | |
quadruple | QReal | |
COMPLEX | single | Complex |
double | DComplex | |
double-double | DDComplex | |
quadruple | QComplex |
DDRealとQRealはほとんどのシステムで同一です。ただしSun Solarisにおいてはこれら は異なります(後者はIEEE互換の128ビット浮動小数点型となります)。
12.1.4.1 ポインタ
スカラで多相でなく文字型でもないPOINTER
型は単にオブジェクトに対するCポインタです。
配列POINTER
型は対象の配列を記述するドープベクトル(dope vectors)です。
単なるアドレスとは異なり、これらのドープベクトルは連続な配列中における不連続
な部分配列を直接記述することができます。
詳細については下記(Dope
N とChDope
N)をご参照ください。
多相のドープベクトルはCSDope
NであるCLASS(*)
ポインタを除き NPDope
N です。
12.1.4.2 構造型
Fortranの構造型はCのstruct
に翻訳されます。
BIND(C)
またはSEQUENCE
が使用された場合には、構造体内の項目の順序は
構造型定義中における順序と同一になります。
そうでない場合には、並び方を改善し記憶容量を減らすべく順序の入れ替えが行われ、
より大きな型が構造体の先頭の方に置かれることになります。順序の確認はCの出力コードを調べる事で行えます。
12.1.4.3 サポートタイプ
Char
unsigned char
基本(単一バイト)文字記憶のデータ型Char2
unsigned short
16ビット(JIS及びUCS-2)文字記憶のデータ型Char4
unsigned int
32ビット(UCS-4)文字記憶のデータ型Substring
struct { Char *addr; Chrlen len;}
単一バイト(基本)CHARACTER
列を記述します。 無指定文字長の基本CHARACTER
変数に使用され、 可変長のスカラnon-POINTER
CHARACTER
関数、 及びすべてのPOINTER CHARACTER
関数の戻りタイプとして使用されます。Substring2
struct { Char2 *addr; Chrlen len;}
16ビット(JISもしくはUCS-2)CHARACTER
文字列を説明し、 無指定文字長の16ビットCHARACTER
変数、 及び可変長スカラ非POINTER
の16ビットCHARACTER
関数の戻り型、 及び全てのPOINTER
16ビットCHARACTER
関数で利用されます。Substring4
struct { Char4 *addr; Chrlen len;}
32ビット ISO_10646(USC-4)CHARACTER
文字列を説明し、 無指定文字長の32ビットCHARACTER
変数、 及び可変長スカラ非POINTER
のUSC-4CHARACTER
関数の戻り型、 及び全てのPOINTER
USC-4CHARACTER
関数で利用されます。Offset
int
,long
またはlong long
アドレス付けと添え字計算のための整数型です。これは32ビットシステム及び小型64ビット システム上ではint
(32ビット)、大型の64ビットシステム上では64ビット整数 型です。Chrlen
- 通常
int
(64ビットのWindowsではlong long
)
文字長を表すための整数型です。 Pointer
- 通常
char *
任意の型を参照するために、及びポインタ算術用に使用されるバイトポインタです。 Triplet
struct { Offset extent,mult,lower;}
配列次元のパラメータを含みます。extent
はその次元のサイズを、mult
は 刻み幅(連続した要素間の距離(単位バイト))を、lower
は下限を表します。 これはDope
N,ChDope
N構造体の成分です。Dope
Nstruct { Pointer addr; Offset offset; Triplet dim[N];}
すべての非多相non-CHARACTER
配列(構造型の配列を含む)に対するドープベク トルです。N は次元数を意味し1から7までです。addr
は最初の要素のアド レスであり、dim
はそれぞれの次元を記述します。offset
は添え字計算 に足されるオフセット、すなわちSUM(mult*lower)
を表します。これはPOINTER
配列関数に対する戻りタイプとして使用されます。それに対するポイン タは形状引継ぎ配列、POINTER
配列引数用の引数タイプとして使用されます。ヌル化された配列ポインタはヌルポインタである
addr
フィールドを持ちます。 ゼロサイズの配列はヌルポインタではないaddr
フィールドを持つので注意し てください。ChDope
Nstruct { Pointer addr; Chrlen len; Offset offset; Triplet dim[N];}
Dope
N 構造体にCHARACTER
長を指定するlen
成分が加わった ものです。ArrayTemp_type
struct { type *addr; Offset extent[7];}
type がInteger1
,Integer2
,Integer
,Integer4
,Logical1
,Logical2
,Logical
,Logical4
,Real
,Double
,QReal
,Complex
,DComplex
,QComplex
のいずれかである連続配列を記述します。 それはnon-POINTER
配列関数に対する戻りタイプとして使用されます。 配列の次元数より後のextent
値は不定である点に注意してください。ArrayTemp_Character
struct { Char *addr; Chrlen len; Offset extent[7];}
連続なCHARACTER
配列を記述します。他の配列型と同じですが、CHARACTER
長に対するlen
成分が追加されています。ArrayTemp_Derived
ArrayTemp_Character
と同義
任意の構造型の連続配列を記述します。 この場合のlen
は構造型配列要素のサイズを表します。
12.1.5 SUBROUTINE戻りタイプ
12.1.5.1 ラベル引数を持つSUBROUTINE
戻りタイプはint
です。その値は制御を移す先のラベルに対するインデックス
(最初のラベルであれば1、以下同様)を表します。
ゼロ、または範囲外の値は制御の移行を伴わないことを示します。
12.1.5.2 ラベル引数を持たないSUBROUTINE
戻りタイプはvoid
です。
12.1.6 FUNCTION戻りタイプ
12.1.6.1
スカラ、non-POINTER
関数
-
INTEGER
,LOGICAL
,REAL
上記組込み型。 -
COMPLEX
, -compatible オプション使用せず(もしくはWindows 64ビットモードではない)
精度に応じてComplex
、またはDComplex
。 -
COMPLEX
, -compatible オプション使用(もしくはWindows 64ビットモード)
戻りタイプはvoid
です。 型がComplex
またはDComplex
の一時データのアドレスが関数の初期引数 として渡されます(結果はその位置に書き込まれます)。 -
CHARACTER
, 固定長または長さ引継ぎ
戻りタイプはvoid
です。 2つの追加の初期引数が関数に引き渡されます。1番目のもの(Char*
,Char2*
,Char4*
のいずれか)は結果 を書き込む場所を示すアドレスを、2番目のもの(Chrlen
)は結果の長さ(呼 ばれた関数が長さ引継ぎ型の場合)を示します。 -
CHARACTER
, 可変長
戻りタイプはSubstring
,Substring2
,Substring4
(上記参照)のいずれかです。 呼ばれた関数は文字列用の記憶領域を割付けます。その領域が必要なくなったときに それを解放するのは呼出し元の責任です。 -
構造型
構造体。
12.1.6.2
スカラPOINTER
関数
すべてのPOINTER
値関数に関し、結果がポイントする記憶領域は呼ばれる関数
内で割付けられていることもあれば、結果が既存の記憶領域をポイントしていること
もある点に注意してください。
-
INTEGER
,LOGICAL
,REAL
,COMPLEX
適切な組込み型(例えばComplex*
)へのポインタ。 -
CHARACTER
戻りタイプはSubstring
,Substring2
,Substring4
のいずれかです。 -
構造型
構造体へのポインタ。
12.1.6.3
配列non-POINTER
関数
-
組込み型
上記組込み型に対する適切なArrayTemp_
構造体。 -
構造型
そのlen
成分が構造体のサイズにセットされた形でArrayTemp_Derived
が 応答として返されます。
12.1.6.4
配列POINTER
関数
-
CLASS(*)
配列の次元によりCSDope1
,CSDope2
, ...CSDope7
のいずれか -
CLASS(
derived-type-name)
配列の次元によりNPDope1
,NPDope2
, ...NPDope7
のいずれか -
非多相の非
CHARACTER
型
配列の次元数に応じてDope1
,Dope2
, ... ,Dope7
のいずれか。 -
CHARACTER
配列の次元数に応じてChDope1
,ChDope2
, ... ,ChDope7
のい ずれか。
非多相の構造型配列はDope
N 構造体として返される点に注意してください。
12.1.7 引数タイプ
12.1.7.1
CHARACTER
型
CHARACTER
型の通常の引数はすべて、CHARACTER
、マルチバイトCHARACTER
を問わず、
CHARACTER
言語要素の長さを示すChrlen
引数を伴っています。この追加引数はすべての通常の引数よりも後、
引数並びの末尾で引き渡されます。CHARACTER
引数が複数ある場合、その長さは左から右に順に引き渡されます。
これについて、32ビットのWindowsの場合で-compatibleオプションを指定した場合を除き、追加の引数はCHARACTER
の直後に続きます。
このルールに対するもう一つの例外は形状引継ぎCHARACTER
配列の場合ですが、この場合、仮引数の長さはドープベクトル中のフィールドから直接取得され、別個に引き渡されることはありません。
12.1.7.2
non-POINTER
非ALLOCATABLE
スカラ
-
non-
CHARACTER
型
引数のアドレスが引き渡されます(例:INTEGER
引数の場合のInteger*
)。 -
CHARACTER
引数のアドレスが引き渡されると共に、引数の長さが別のChrlen
引数として引数並びの末尾で引き渡されます。
12.1.7.3
POINTER
及びALLOCATABLE
スカラ
これらはポインタのアドレス(もしくは割付記述子)が引き渡される点を除くと、通常の場合と同様の形で引き渡されます(例:INTEGER POINTER
の場合のInteger**
)。
12.1.7.4
non-POINTER
配列
-
引継ぎ形状
配列を記述する適切なドープベクトルのアドレスが引き渡されます (すなわち多相及び仮引数の型により、CSDope
N*
,NPDope
N*
,Dope
N*
,ChDope
N*
のいずれか)。 配列は部分配列で記述できる限り連続である必要はありません。 -
その他
配列の最初の要素のアドレスが引き渡されます。CHARACTER
配列の場合、それぞれの配列要素の長さはスカラに対する場合と 同様、引数並びの末尾でChrlen
として引き渡されます。 配列は連続でなくてはなりません。
12.1.7.5
POINTER
配列
引継ぎ形状配列の場合と同様、適切なドープベクトルのアドレスが引き渡されます。
CHARACTER
配列の場合、長さは別個の引数として末尾で引き渡されます。
12.1.7.6 手続き
手続きのアドレスが引き渡されます。CHARACTER
関数の場合、長さは別個の引数として並びの末尾で引き渡されます。
関数が可変長型の場合、この長さは負の値となります。
12.1.7.7
OPTIONAL
引数
OPTIONAL
な引数が.NOT.PRESENT()
である場合、適切な型のヌルポイン
タが引き渡されます(例えばINTEGER
スカラの場合、(Integer*)0
が
引き渡されます)。